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ABSTRACT 


The  motivation  for  this  study  lies  in  the  significance  of  free  surface  suction 
effects  during  submarine  operations  at  periscope  depth.  Such  operations  become 
increasingly  important  as  new  roles  for  the  Navy  in  littoral  waters  are  emerging. 
Particular  emphasis  is  placed  on  computation  of  steady  state  forces  on  the  body  as 
a  function  of  speed,  depth,  and  wave  frequency  and  direction.  These  forces 
constitute  an  important  and  very  frequently  limiting  factor  in  establishing  the 
periscope  depth  submerged  operating  envelope.  Solution  of  the  problem  is 
accomplished  by  singularity  distribution  on  the  actual  surface  of  the  body  and 
discretization  in  the  form  of  plane  quadrilateral  elements.  Parametric  studies  are 
conducted  in  order  to  assess  the  effects  of  body  shape  and  size.  The  results  of  this 
thesis  can  be  directly  utilized  in  the  simulation  based  design  process  as  well  as 
during  training.  [ 
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I.  INTRODUCTION 


S\ibinarine  missions  are  rapidly  evolving  into  frequent 
deployments  to  littoral  areas,  which  require  operating  in 
shallow  water.  Shallow  water  is  usually  defined  as  a  charted 
depth  of  less  than  600  feet,  but  can  be  as  shallow  as  180  feet 
or  less  in  certain  situations.  Shallow  water  operations 
sharply  contrast  with  the  traditional,  blue  water  missions  of 
the  Cold  War.  Instead  of  concentrating  on  anti-submarine 
warfare  (ASW)  against  modem  Soviet  platforms,  submarines  will 
execute  ASW  against  older,  diesel  submarines  operating  in 
littoral  areas.  Also,  surveillance  (both  visual  and 
electronic  )  operations,  mine  laying  evolutions,  and  joint 
special  forces*  escpeditions  will  become  routine  occurrences  as 
the  submarine  force  responds  to  dynamic  regional  conflicts. 

Given  the  edsundant  time  submarines  will  spend  operating  in 
shallow  water,  the  importance  of  periscope  depth  operations 
greatly  expands.  Periscope  depth  is  the  depth  at  which  a 
submarine  interacts  with  the  outside  world.  It  has  a 
periscope  exposed  along  with  one  or  more  multi-purpose  masts. 
During  routine  operations,  a  submarine  is  at  periscope  depth 
to  conduct  communications  (receive  and  transmit) ,  to  establish 
navigational  fixes,  to  perform  housekeeping,  and  to  ventilate 
the  submarine  atmosphere.  During  operations  in  shallow 
waters,  the  time  spent  at  periscope  depth  will  greatly 


1 


increase  and  may  be  the  only  depth  the  sxibmarine  is  edsle  to 
operate.  As  a  result,  the  hydrodynamic  influences  on  the 
submarine  at  periscope  depth  are  of  special  interest. 

This  work  primarily  examines  the  free  surface  effects  on 
submerged  bodies  similar  to  submarine  hulls.  As  a  submerged 
body  travels  through  the  water,  it  generates  waves  which 
interact  with  the  free  surface.  This,  in  turn,  produces 
forces  and  moments  on  the  body.  A  longitudinal  force,  drag, 
opposes  the  body's  motion,  while  a  normal  force,  lift,  pushes 
the  body  to  the  surface.  Additionally,  a  moment  is  produced 
on  the  body  which  causes  the  body  either  to  pitch  up  or  down. 
Of  these  three  effects,  lift  is  the  most  important,  since  it's 
desircible  to  keep  the  body  from  breaking  the  free  surface 
(broaching) .  Broaching  hazards  the  submarine  by  making  it 
much  more  susceptible  to  detection  from  air,  sea  or  land 
sources.  It  is  also  desiredsle  to  minimize  the  moment  in  order 
that  the  stem  (propeller)  not  break  the  free  surface.  A 
propeller  out  of  the  water  greatly  reduces  the  submarine's 
speed,  making  depth  control  much  more  difficult  and  also 
increases  the  sxibmarine's  detectcdsility.  Prior  to  proceeding 
to  periscope  depth,  a  submarine  crew  will  adjust  its  trim  by 
taking  on  additional  weight  (seawater)  into  forward,  aft  or 
auxiliary  (middle)  ballast  tanks.  This  ballast  will  help 
counteract  the  upwards,  sea- suction  lift  force  encountered 
when  operating  near  the  free  surface.  However,  taking  on  too 
much  ballast  makes  it  difficult  to  reach  periscope  depth  and 
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could  cause  the  sxjbmarine  to  lose  depth  control  (sink  out) 
once  periscope  depth  is  reached.  Once  at  periscope  depth,  the 
ordered  depth  is  maintained  using  a  combination  of  speed, 
control  surfaces  (bow/fairwater  planes,  stem  planes,  and 
rudder)  and  variable  ballast.  These  are  used  to  prevent 
broaching  and  to  keep  the  propeller  submerged. 

In  addition  to  the  waves  the  submerged  body  produces,  the 
existing  seas  greatly  effect  the  operation  of  the  body  near 
the  surface.  A  submarine  can  easily  maintain  depth  in  calm 
seas,  but  in  the  presence  of  high  sea  states,  depth  control  is 
much  more  difficult.  The  submarine  pitches  and  rolls  in  the 
presence  of  sea  waves  and  this  reduces  the  crew's  edsility  to 
keep  the  submarine  submerged.  The  resulting  drag,  lift  and 
moment  from  the  incident  waves  are  presented  in  this  work. 
The  results  from  the  wave  making  and  the  incident  wave 
problems  can  be  linearly  combined  to  produce  a  total  effect  on 
the  submerged  body. 

On  the  practical  side,  the  results  of  this  work  can  be 
utilized  in  the  simulation  based  design  process  as  well  as 
during  training  in  sxibmarine  simulators.  Additionally,  they 
could  be  incorporated  into  future  automatic  depth  control 
system  designs.  Also,  incorporation  into  the  existing  SFMPL 
(Ship's  Fleet  Mission  Program  Library)  would  allow  a  submarine 
crew  to  make  accurate  predictions  of  the  ballast  required  at 
periscope  depth.  The  wave  making  and  incident  wave  results 
are  linearly  added  to  the  forces  and  moments  produced  by  other 
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significant  effects  such  as  added  znass  from  the  body's 
acceleration,  frictional  resistance  due  to  viscosity,  and  eddy 
resistance  from  appendages,  to  create  a  combined  total  force 
and  moment  acting  on  the  body.  The  addition  of  these  two 
effects  due  to  operation  near  the  free  surface  will  greatly 
enhance  the  accuracy  and  realism  of  the  simulators  and  iit^rove 
future  sxibmarlne  designs. 

For  this  work,  the  Defense  Advcuiced  Research  Projects 
Agency  (DARPA)  SUBOFF  model  (DTRC  Model  5470)  (Roddy, 1990)  is 
utilized  to  provide  the  sizbmerged  hull  shape  into  the  computer 
program.  The  FORTRAN  program  written  by  Doctors  and  Beck 
calculates  the  drag,  lift  and  moment  coefficients  for  a  given 
set  of  input  conditions.  For  the  wave  making  or  Neumann- 
Kelvin  problem,  the  free  surface  flow  created  by  a  moving 
submerged  body  is  solved  for  (Doctors  &  Beck, 1987).  The 
program  yields  a  solution  to  this  problem,  in  which  the  fluid 
is  inviscid  and  exact  body  boundary  conditions  (no  normal 
velocity  on  the  body  surface)  are  satisfied.  The  sea  surface 
condition  is  linearized,  rather  than  being  met  exactly.  The 
SUBOFF  model  is  discretized  into  panels  of  constant  source 
strengths  using  the  method  of  Hess  and  Smith  (Hess  & 
Smith, 1964  and  Parsons, 1984) .  The  source  strengths  are  then 
solved  for  such  that  the  body  boundary  condition  is  satisfied. 
For  the  incident  wave  problem,  the  progreun  uses  a  given 
es^resslon  for  the  pressure  distribution,  applies  it  to  each 
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discrete  hull  panel  and  then  svuns  the  resulting  forces  euid 
moment.  This  is  examined  for  various  depths,  body  speeds  and 
vrave  directions. 

In  this  work,  the  two  problems  are  initially  formulated  in 
Chapter  II.  The  Inviscid  flow  theory  development  and  cui 
explanation  of  the  boundary  conditions  es^lain  the  theoretical 
background  needed  to  understand  the  wave  making  problem 
(Doctors  &  Beck, 1987  and  Papoulias, 1993  and  Papoulias  & 
Beck,  1988  and  Reed  &  Beck  &  Griffin  &  Peltzer,  1990)  .  The 
incident  wave  problem  is  next  set  forth  (Papoulias, 1993) 
followed  by  a  brief  discussion  on  the  niimerical  solution 
techniques.  Chapter  III  presents  the  results.  After 
determining  the  number  of  integration  points  needed  for 
convergence,  the  results  of  the  wave  making  and  incident  wave 
problems  are  presented.  The  wave  making  runs  show  the  drag, 
lift  and  moment  as  a  function  of  both  speed  and  depth.  The 
incident  wave  runs  show  these  three  parcuneters  as  a  function 
of  depth,  body  speed  and  wave  direction  for  various 
wavelengths.  After  exploring  the  wave  making  and  incident  wave 
problems.  Chapter  IV  presents  a  parametric  study  of  submarine 
hull  shapes.  The  siUsmarine  hull  is  described  by  three 
sections,  a  cylindrical  parallel  midbody  with  a  parcJsola  of 
revolution  bow  (entrance)  and  an  ellipsoid  of  revolution  stem 
(run)  (Jackson, 1992) .  The  shape  of  the  entrance  and  run  can 
be  varied  by  changing  coefficients  in  the  mathematical 
expression  for  each.  Two  effects  are  examined.  In  the  first. 
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the  displacement  cuid  the  diameter  are  fixed  and  the  length  is 
varied.  In  the  second  case,  the  displacement  and  the  length 
remain  constant  and  the  diameter  is  changed.  In  both  cases, 
the  drag,  lift  and  moment  are  investigated  for  both  the  wave 
making  and  incident  wave  cases.  These  parametric  studies  have 
direct  application  in  the  design  process  of  future  classes  of 
submarines.  For  a  given  set  of  rec[uirements,  the  effect  on 
lift  and  moment  can  be  e^qplored  to  provide  an  optimum  solution 
to  periscope  depth  operations.  Thus,  this  work  increases  the 
understanding  of  near  surface  effects  on  a  submerged 
submarine . 
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II.  PROBLEM  FORMULATION 

A.  INTRODUCTION 

In  the  development  of  the  free  surface  effects  on 
submerged  bodies,  the  forces  and  moments  acting  on  the 
submerged  body  must  calculated.  The  forces  include  the  drag 
(force  parallel  to  the  body) ,  the  lift  (force  norxnal  to  the 
body)  and  the  bow-up  moment.  These  hydrodynamic  forces  acting 
on  the  submerged  body  arise  from  a  modification  to  the 
pressure  distribution  summed  around  the  surface  area  of  the 
body  in  question.  Forces  can  arise  from  relative  motion 
velocity  between  the  body  and  the  fluid.  This  relative  motion 
can  be  either  time  varying  or  time  invariant.  The  approach 
utilized  here  considers  only  the  time  invariant  cases  of  the 
motion  of  the  body  through  a  stationary  fluid  or  the  motion  of 
the  seaway  relative  to  the  body.  Thus,  the  inertia  forces 
from  an  accelerating  body  or  from  time  varying  seaways  are 
neglected.  This  assumption  siirplifies  the  problem  and  is  a 
first  step  towards  a  comprehensive  modelling  of  a  s\ibmarine  at 
periscope  depth. 

In  examining  the  combined  effects  of  the  relative  motion 
between  the  body  euid  the  fluid,  the  system  is  considered  to  be 
a  linear  system  and  thus  the  principle  of  superposition 
applies.  Separate  effects  can  be  examined  independently  and 


the  results  summed  together.  This  work  considers  two  separate 
effects.  The  first,  the  wave  incdclng  problem,  Is  the  motion  of 
the  body  relative  to  a  stationary  fluid,  resulting  in  forces 
and  moments  being  generated  on  the  body.  The  free  surface  is 
assumed  to  be  an  otherwise  calm  sea.  Any  modification  in  the 
free  surface  shape  (waves)  is  due  to  the  influence  of  the  body 
itself.  The  second  consideration  is  the  effects  of  an 
existing  seaway  on  the  body  or  the  incident  wave  problem. 
These  two  effects  can  then  be  summed  together  to  yield  the 
total  response  of  both  effects. 

The  appropriate  fluid  mechanics  tools  that  are  used  to 
describe  sea  waves  and  ship  motions  are  based  on  potential 
flow  theory.  The  first  fundamental  assuirqption  of  ideal 
(potential)  flow  theory  is  that  mass  is  preserved.  If  we 
consider  a  control  volume  surrounding  the  fluid  and  bodies  of 
interest,  what  mass  enters  the  volume  either  accumulates 
inside  or  leaves  the  control  vol\une.  This  can  be  esqpressed 
mathematically  by  utilizing  the  divergence  theorem:  Given  a 
closed  surface  area  S  with  a  unit  vector  h  pointing  inward, 
and  volume  V  enclosed  by  S,  then  for  any  single  valued  and 
differentiable  vector  function  A, 

JJ^JT-UdS  «  -  JJJ^V-AdV  .  (1) 
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If  we  denote  the  xanit  vector, 


B  -  +  ziyS  +  njc 


(2) 


then. 


A  •  fi  =  AjjDjp  +  AyZly  +  A^n,  , 


(3) 


which  is  the  component  of  A  in  the  direction  of  !!.  If  we 
denote  the  operator, 

V  =  JLI  +  4-5  +  J-k  ,  (4) 

Tx  Tx-'  Tz 


then  the  dot  product. 


_  _  dAx  9Ay  BA, 

V*A  =  -— ir  +  -  y  + 

Bx  By  Bz 


(5) 


Conservation  of  mass  then  requires. 


(6) 


where  the  right  hand  side  represents  the  net  mass  flow  in,  and 
the  left  hand  side  the  increase  of  mass  in  V.  Applying  the 
divergence  theorem,  we  get, 

,  (7) 
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and,  since  V  is  an  arbitrary  volume. 


bp 

Ti 


+  V-  (pff) 


0 


(8) 


If  we  assiime  that  the  density  p  is  consteuit,  a  good 
assumption  in  naval  hydrodynamics,  we  get, 

V •  U  ‘  0  . 


If  the  velocity  vector  U  is, 

fj  -  ul  +  *  wfc 


(10) 


we  get  the  final  form  of  the  continuity  equation. 


bu 

T5t 


+ 


i 


(11) 


inside  the  control  volume  V. 

The  next  very  important,  and  not  so  obvious,  assurr^jtion  is 
that  the  flow  is  irrotational  (flow  without  rotation) .  One 
important  property  of  irrotational  flows  is  that  the 
circulation  around  any  closed  curve  c  is  zero, 

dS  •  0  .  (12) 


This  means  that  this  line  integral  will  be  independent  of  the 
path  of  integration,  which  can  only  be  achieved  if  dS’U  is  an 
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we  get, 

or,  in  vector  form, 

U  ^V<f>  ,  (18) 

which  is  the  definition  of  the  velocity  potential.  This  is  a 
very  useful  quantity  since  instead  of  confuting  the  fluid 
velocity  which  is  a  vector  function,  0,  all  we  have  to  do  is 
compute  the  scalar  velocity  potential  0,  and  then  we  can  get 
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the  components  of  0  by  differentiation.  Using  4>i  we  can  write 
the  continuity  equation  as  follows, 


or. 


^JCC  *  4>yy  *  4>ZZ  *  0 


(20) 


which  is  Laplace's  equation. 

This  is  the  equation  we  have  to  solve,  subject  to 
appropriate  boundary  conditions,  to  compute  the  flow  field  of 
an  ideal  flow.  Once  we  compute  e^nd  the  velocities,  we  can 
find  pressure  by  using  Bernoulli's  ecjuation, 

p  +  +  pgz  *  c  ,  (21) 

where  c  is  a  constant,  taken  equal  to  zero  in  most  cases. 
This  form  is  similar  to  the  usual  form  of  Bernoulli's 
equation,  with  the  extra  term  d4>/dt  due  to  the 

possible  unsteady  (always  irrotatlonal)  nature  of  the  flow. 
The  pressure  p  can  then  be  computed  from, 

p  »  -p-||  -  lpV4>-V4>  -  pgz  .  (22) 


12 


The  first  two  terms  represent  the  hydrodynamic  contribution  to 
the  pressure  and  the  third  term  represents  the  hydrostatic 
contribution.  Integration  of  the  pressure  distribution  then 
over  the  surface  of  the  body  will  produce  forces  and  moments. 

To  summarize,  the  basic  features  of  potential  flows  are: 

1.  Mass  is  conserved  eind  fluid  density  is  constant.  This 
results  in, 

V-  a  »  0  ,  ^23) 


the  continuity  equation. 

2.  The  flow  is  irrotational ,  which  means  we  can  define  the 
velocity  potential, 

U  =  ,  (24) 

which  satisfies  Laplace's  equation, 

^XX  *  4>yy  *  4>ZZ  *  ^  *  (25) 

3.  For  unsteady,  irrotational  flows,  we  can  use  Bernoulli's 
equation, 

p  +  +  ipV^*V0  +  pgz  =  0  .  (26) 

This  comes  from  Newton' s  law,  and  as  such  it  is  relative  to  an 
inertial  coordinate  system.  It  is  noted  here  that  throughout 
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the  remainder  of  this  work,  partial  derivatives  will  be 
designated  by  a  subscript  vice  the  d/dx  notation  (i.e., 
vice  d^/dx) . 

Wave  effects  on  the  free  surface  are  iir5)ortant  effects  to 
consider  on  a  submerged  body.  The  most  common  cind  in^ortant 
wave  phenomena  are  free  surface  waves  that  exhibit  typical 
periods  of  a  few  seconds.  Other  waves  like  sxibsurface  or 
internal  waves  are  found  in  regions  of  density  stratification 
beneath  the  free  surface  and  are  typically  of  lower  frequency 
with  periods  on  the  order  of  several  minutes.  The  influence 
of  such  waves  on  submerged  bodies  is  generally  negligible 
unless  the  body  has  an  unusually  low  frequency  resonance.  The 
simplest  free  surface  wave  formation  is  the  plane  progressive 
wave  system.  This  motion  is  two  dimensional,  sinusoidal  in 
time  with  angular  frequency,  w,  and  propagates  with  phase 
velocity,  Cp,  such  that  to  an  observer  moving  with  this 
velocity  the  wave  appears  to  be  stationary.  The  waves  could 
also  be  modelled  as  a  series  of  random  waves,  where  random 
refers  to  the  character  of  the  wave  height  distribution. 
These  can  be  represented  using  a  probedsilistic  approach.  This 
work,  however,  views  the  waves  as  plane  progressive  waves  of 
small  amplitude,  with  sinusoidal  time  dependence.  Effects  of 
changing  depth,  varying  wave  period  (and  frequency)  and 
altering  wave  direction  are  explored.  Additionally,  the 
interaction  of  the  incident  waves  on  the  body,  known  as  wave 
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diffraction,  is  neglected.  This  is  a  valid  assumption  since 
the  incident  wavelengths  are  the  same  magnitude  as  the  body 
length  and  significant  diffraction  is  not  expected  to  occur. 


B.  WAVE  MAKING  PROBLEM 

In  the  following,  a  right  handed  coordinate  system  is  used 
with  the  origin  at  the  projection  onto  the  free -surface  of  the 
intersection  of  the  horizontal  centerline  and  midships.  The 
positive  x-aocis  points  out  the  bow,  the  y-axis  is  positive  to 
port  cind  the  z-axis  in  position  upward.  The  submerged  body  is 
at  a  depth  H  below  the  undisturbed  free  surface  and  of 
diameter  D  at  midships.  Figure  1  illustrates  this. 


The  total  potential  for  a  sxibmerged  body  advancing  with  a 
constant  speed  U  in  an  otherwise  calm,  inviscid  and 
incompressible  fluid  is  given  by 

♦  (x,y,  z)  •  -  Ux  *  i>{x,y,z)  ,  <27) 

where  ♦  is  the  free  stream  potential  and  ^  is  the  perturbation 
potential  due  to  the  motion  of  the  body  alone.  Assuming  is 
small  relative  to  Ux,  the  boundary  value  problem  and  its 
boundary  condition  may  be  rewritten  in  terms  of  As  such, 
<t>  satisfies  the  Laplace's  equation 

V2<|>  =  0  (28) 

On  the  free  surface,  two  boundary  conditions  apply.  The 
first  is  the  dyneimic  free -surface  condition,  which  ensures 
constant  pressure  on  the  free  surface.  The  wave  elevation, 
i^(x,y),  may  be  expanded  in  a  Taylor  series  about  z-O,  where 
all  higher  order  terms  are  neglected  to  obtain 

Sr7  +  (V^  '  on  z  -  i»(x,y) 

where  g  is  the  acceleration  due  to  gravity.  The  second  is  the 
kinematic  free-surface  boundary  condition,  which  requires  that 
there  be  no  flow  through  the  free-surface. 


*xVx  *  ♦yV  “  ♦*  “  0  ,  on  z  =  ti{x,y)  .  (30) 

Rewriting  the  dynamic,  (29) ,  and  the  kinematic,  (30) ,  free- 
surface  conditions: 

Sn?  +  ^{2U<l>jc  +  (^jc)^  +  +  (^*)^}  *  0  on  z  «  0  ,(31) 

and 

^Vx  *  ^xVx  "*■  4>yfly  ”  ♦z  “  0  on  z=0  .  (32) 

Neglecting  the  quadratic  terms  in  both  equations,  and 
differentiating  the  dynamic  free- surface  condition  with 
respect  to  x  and  subtracting  the  result  from  the  kinematic 
conditions,  gives  the  linearized  free-surface  condition: 

4>xx  *  “0  ^“0  ,  (33) 


where  ko  ■  g/U^. 

The  boundary  condition  on  the  surface  of  the  body,  the 
body  boundary  condition,  states  that  there  shall  be  no  flow 
through  the  surface  of  the  body: 

^  *  0  ,  on  f  (x,y)  -  z  =  0  ,  (34) 

where  d/dn  denotes  the  derivative  in  the  direction  of  the 
three  dimensional  normal  vector  pointing  into  the  body  (h)  and 
f  (x,y)  represents  the  surface  of  the  body.  This  body  boundary 
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condition  can  also  be  written  as 


d<l> 

IB 


s  •  -  u. 


/ 


1  +  +  (fv) 


on  f{x,y)  -  z  *  0 

(35) 


A  fundamental  solution  of  (28)  and  (33)  is  the  Green 
function,  given  by 

GiP.Q)  =  1  -  J,.  +  S{P,Q)  ,  (36) 

r 


where 


GiP,Q) 


^oLim  rw  (»^j^exp{k[z+t-^i{x-i)cos6*i  (y-T?)sing] } 
TT  Jo  _  jccos^®  -  i/iCOS® 


(37) 


where  P  ■  (x,y,z)  is  the  field  point,  Q  ■  (iiViH)  is  a  point 
source  of  strength  -4ir  on  the  body  ,  and 

r,  r'  «  [(x-{)^  +  (y-1))^  +  (zT 

In  this  definition,  r  is  the  distance  between  the  field  point 
and  the  singularity  point  and  r'  is  the  distance  between  the 
field  point  and  the  image  of  the  singularity  point  in  the 
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free -surface.  The  denomiziator  of  the  integrand  in  (37) 
contains  the  Rayleigh  virtual  viscosity  fi,  which  is  taken  to 
be  sxnall  and  positive,  thus  ensuring  that  the  radiation 
condition  is  satisfied.  The  1/r  terms  in  (36)  represent 
source  and  sink  distributions  typical  in  potential  flow 
problems.  The  extra  term,  (37) ,  represents  a  series  of  waves 
such  that  the  free  surface  boundary  condition  (33)  is 
automatically  satisfied,  where  6  is  the  wave  direction  emd  k 
is  the  wavenumber.  Finally,  the  radiation  condition  (i.e., 
all  waves  far  away  from  the  body  are  outgoing)  is  satisfied  by 
taking  /x  to  be  positive.  The  only  condition  that  still  needs 
to  be  satisfied  is  the  body  boundary  condition  (35) . 

In  the  usual  method  of  potential  theory.  Green's  theorem 
can  be  applied  to  a  large  voltune  of  fluid  containing  the  body, 
and  extending  to  infinity  both  laterally  and  in  depth.  The 
following  result  for  the  perturbation  potential  in  terms  of  a 
pure  source  strength,  a,  is  obtained: 

<t>{P)  -^jj^G(P,Q)o{Q)dS{Q)  ,  (39) 


where  Su  is  the  wetted  surface  of  the  body. 

An  integral  equation  for  the  source  strength  may  be  found 
by  differentiating  (39)  with  respect  to  the  normal  on  the  body 
and  setting  it  equal  to  the  body  boundary  condition  which 
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requires  that  the  normal  velocity  on  the  hull  be  zero. 

Using  (27)  and  (39),  we  may  write 

-^a-^[\^Ga(P,Q)a{Q)<iS{Q)  .  (40) 

Solution  of  (40)  is  achieved  by  discretizing  the  surface 
of  the  body  into  plane  quadrilateral  panels  using  the  method 
of  Hess  and  Smith,  as  described  later.  The  equations  are  then 
assembled  to  yield  a  system  of  linear  algebraic  equations  in 
terms  of  the  unknown  source  a.  This  method  of  solution  is 
also  known  as  the  Neumann- Kelvin  method.  The  flow  past  a 
body  moving  at  a  steady  speed  requires  the  body  boundary 
condition  be  satisfied  exactly,  while  the  free-surface 
condition  is  satisfied  in  a  linearized  sense. 

C.  INCIDENT  NAVES  PROBLEM 

A  second  major  effect  on  a  submerged  body  operating  near 
the  free  surface  is  that  of  incident  waves  disturbing  the 
otherwise  calm  free  surface  and  interacting  with  the  body.  In 
the  case  examined  in  this  work,  the  waves  incident  upon  the 
body  are  two  dimensional  plane  progressive  waves  of  small 
eunplitude,  with  sinusoidal  time  dependence.  The  wave  motion 
is  parallel  to  the  x-z  plane,  with  angular  frequency  «, 
propagating  with  phase  velocity  Cp.  The  body  motions  are 
assumed  to  be  sufficiently  small  so  that  linear  theory  holds. 
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The  free  surface  elevation  Is  of  the  general  form 


ij  (x,  t)  *  Acoaikx  -ut) 


(41) 


where  the  positive  x-axis  is  chosen  to  coincide  with  the 
direction  of  wave  propagation.  Here  A  is  the  wave  amplitude, 
and  the  parameter  k  «  w/Cp  ,  is  the  wavenumber,  the  number  of 
waves  per  unit  distance  along  the  x-axis.  Clearly  k  -  2ir/X  , 
where  the  wavelength  X  is  the  distance  between  successive 
points  on  the  wave  with  the  same  phase. 

The  solution  of  this  problem  is  expressed  in  terms  of  a 
two  dimensional  velocity  potential  ^(x,z,t)  which  must  satisfy 
Laplace's  equation 

V20  =  0  ,  (42) 


and  appropriate  boundary  conditions.  Furthermore,  ^  must 
yield  the  wave  elevation  (41)  from 


(43) 


Equation  (43)  is  the  so  called  linearized  dynamic  boundary 
condition  on  the  free  surface  and  is  an  expression  of  the 
fact,  through  Bernoulli's  equation,  that  the  pressure  on  the 
free  surface  must  be  the  same  as  the  ambient  atmospheric 
pressure. 
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An  appropriate  boundary  condition  on  the  sea  bottom  is 

*  0  ,  at  z  «  -h  ,  (44) 

i.e.,  the  bottom  at  depth  h  is  a  rigid  ixrpermeable  plane. 
Finally,  the  free  surface  boundary  condition  is 

*  9^z  *  0  »  on  2*0  .  (45) 

Equation  (45)  is  a  combined  dynamic  and  kinematic  free  surface 
boundary  condition.  The  dynamic  condition  was  mentioned 
before,  while  the  kinematic  condition  simply  states 

Vt  * 

i.e.,  the  vertical  velocities  of  the  free  surface  and  fluid 
particles  are  the  same.  Combining  (43)  and  (46)  we  arrive  at 
(45) ,  ignoring  the  small  departures  of  the  free  surface  from 
the  horizontal  orientation  z  -  0. 

Clearly  the  velocity  potential  4>  must  be  sinusoidal  in  the 
same  sense  as  (41) ;  therefore  we  seek  a  solution  of  the  form 

4>{x,z,t)  -9e{z(z)e-^^*^“‘=}  .  (47) 

Substituting  (47)  into  (42),  Z  must  satisfy  the  ordinary 
differential  equation 


22 


Z„  -  k=Z  .  0  .  (48) 

throughout  the  domain  of  the  fluid.  The  most  general  solution 
of  (48)  is  given  in  terms  of  esqponential  fiinctions  in  the  form 

Z  ■  Ce**  +  Pe”**  .  (49) 

Assuming  infinite  water  depth  (deep  water) ,  the  constant  D  in 
(49)  must  be  zero  to  avoid  an  unbounded  motion  deep  beneath 
the  free  surface,  resulting  in 

Z  =  Ce**  ,  (50) 

Substituting  into  (47) 

^  =  3e{Ce**‘^***^"*=}  .  (51) 

Now  substituting  (51)  into  (43)  with  z  «  0,  and  comparing  the 
result  with  (41)  we  can  find 

C  =  iMh  ,  (52) 

0) 


and 

4>  “  J^e**sin(iac  -  «t) 

Cl) 


(53) 
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An  additional  relation  between  the  wavenumber  k  and  the 
frequency  «  can  be  obtained  by  substituting  (53)  into  (45) . 
This  relation,  called  a  dispersion  relation  is 


k  = 


w' 


9 


(54) 


The  frequency  w  can  be  replaced  by  the  wave  period  T  -  27r/a), 
just  as  k  can  be  replaced  by  the  wavelength  X  «  2ir/k. 

The  phase  velocity  Cp  can  be  determined  from  the 
definition  of  the  wavenumber  and  (54)  to  yield 


Cl) 


(55) 


Equation  (55)  states  that  surface  waves  in  deep  water  are 
dispersive;  longer  waves  travel  faster  than  shorter  waves. 

While  the  wave  moves  with  the  phase  velocity  Cp,  the  fluid 
itself  moves  with  a  much  smaller  velocity  given  by  the 
gradient  of  the  potential  (53)  .  The  velocity  components  (u,w) 
of  the  fluid  are 

u  =  =  wAe**cos  {kx  -  ut)  ,  (56) 


w  ^  (f>^  ^  ajAe**sin(kx  -  ut) 


(57) 
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within  linear  theory,  the  fluid  particles  move  in  small 
circular  orbits  proportional  to  the  wave  amplitude;  they 
remain  in  the  same  mean  position  as  the  wave  propagates 
through  the  fluid  with  a  phase  velocity  independent  of  the 
wave  arplitude.  As  the  depth  of  submergence  beneath  the  free 
surface  increases,  the  fluid  velocities  (56)  euad  (57)  are 
attenuated  exponentially.  For  a  submergence  of  half  a 
wavelength,  kz  >  -x,  the  exponential  factor  is  reduced  to 
0.04.  Thus,  waves  in  deep  water  are  confined  to  a  relatively 
shallow  layer  near  the  free  surface,  with  negligible  motion 
beneath  a  depth  of  ed^out  0.5X.  On  this  basis  one  can 
anticipate  that  if  the  fluid  depth  is  finite,  but  greater  than 
half  a  wavelength,  the  effects  of  the  bottom  will  be 
negligible. 

For  a  fluid  at  constant  depth  h,  the  boundary  condition 
(44)  is  imposed.  Returning  to  the  general  solution  (49) ,  both 
exponential  functions  are  retained,  with  the  constants  C  and 
D  suitably  chosen  to  satisfy  (44)  .  In  this  case,  the  velocity 
potential  (53)  becomes 

<l>  =  i[^£?gMf^^.,:!llsin(kx  -  at)  ,  (58) 

b)  coshkh 

and  the  dispersion  relation  (54)  is  now 

k  tanh  (kh)  «  ~  ,  (59) 

9 
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where  we  Ceui  see  that  as  h  -»  «  we  get  the  previous  deep  water 
results.  The  fluid  velocity  con^onents  cam  be  coii^)uted  as  in 
(56)  eind  (57) ,  and  for  finite  depth  it  follows  that 

u  -  ^  *...i-2eos  (Jbc  -  «t)  ,  (60) 

(I)  cosnxfi 

w  .  -  «t)  .  (61) 

w  coshxfi 

In  this  case  the  fluid  particle  trajectories  are  elliptical, 
with  the  xnajor  axis  being  horizontal. 

The  phase  velocity  for  finite  depth  can  be  expressed  in 
the  form 

Cp  =  ^  =  J  ^tantikh  . 

This  tends  to  the  deep  water  limit  (55)  for  kh  >  1.  The 
opposite  limit,  where  kh  <  1,  is  the  regime  of  shallow  water 
waves.  In  this  case  (62)  can  be  approximated  using  the  Taylor 
series  expansion  for  the  hyperbolic  tangent,  and  the  leading 
order  approximation  for  the  phase  velocity  is 

Cp  *  ,  for  kh<  1  . 
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In  the  shallow  limit,  the  resulting  wave  motion  is 
nondispersive,  since  the  phase  velocity  depends  only  on  the 
depth. 

Once  the  velocity  potential  is  determined  we  can  obtain 
the  total  pressure  by  substitution  into  Bernoulli's  Equation 
(22)  which  is  repeated  here 

P  “  ~P4>t  ” 


Using  the  deep  water  expression  for  the  velocity  potential  0 
(53)  and  solving  (64)  results  in 

p  «  -pgAe**cos  (Jcx  -  wt)  -  -  pgz  .  (65) 

The  first  term  in  (65)  is  the  time  rate  of  change  of  the 
velocity  potential  and  is  neglected  in  this  work  since  only 
time  invariant  effects  are  considered.  The  second  term  in 
(65)  is  the  key  element  examined  in  this  work.  It  represents 
the  time  invariant  effect  on  the  body  and  is  the  major 
contributor  to  the  sea- suction  forces  on  the  body.  The  last 
term  in  (65)  is  the  hydrostatic  coitponent  of  the  total 
pressure  and  is  likewise  neglected.  As  before,  the  pressure 
is  integrated  over  the  entire  surface  to  yield  the  resulting 
forces  and  moments. 
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The  previous  equations  of  motion  are  valid  for  a  body  with 
zero  forward  speed.  If  the  body  possesses  a  forward  speed  U, 
this  can  be  assumed,  within  linearity,  constant.  The  only 
change  in  such  a  case  is  in  the  frequency  w  due  to  a  Doppler 
shift  effect.  If  we  consider  the  case  of  head  seas,  or  waves 
from  directly  ahead,  we  can  assume  that  both  the  wave 
excitation  forces  and  the  resultant  oscillatory  motions  are 
linear  and  harmonic,  acting  at  the  frequency  of  wave  encounter 

2 

-  CD  +  «  CD  +  JtD  -  CD  +  —U  .  (66) 

*  Tfa  gr 

To  account  for  the  more  general  cases  of  waves  at  a  direction 
B,  where  B  >  180^  corresponds  to  head  seas  and  0  >  0**  to 
following  seas,  the  frequency  of  encounter  cd^  becomes 

2 

CD-  =  CD  -  —UcobB  *  CD  -  kUcoaB  .  (67) 

9 

The  frequency  of  encounter  and  wave  direction  are  utilized  in 
this  work  to  determine  the  incident  waves  effect  on  the 
submerged  body.  It  should  be  enphasized  that  the  mean  second- 
order  forces  acting  on  a  near  surface  submarine  in  a  seaway 
are  due  entirely  to  time-averaged  products  of  first-order 
quantities  at  each  separate  frequency  of  monochromatic  wave 
encounter.  In  general  these  mecui  forces  involve  quadratic 
products  of  the  radiation  and  diffraction  potentials.  The 
former  is  usually  neglected  (Bingham  &  Korsmeyer  & 
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Newman, 1994) ,  which  corresponds  to  the  case  where  the 
submarine's  unsteady  motions  are  neglected. 

Thus,  the  mean  second  order  force  can  be  written  as 

F(«)  -  -^pjj^dSnV{(t>j  *  4>s) 

Here  the  total  diffraction  potential  (0j  +  ^g)  is  for  a 
steady- state  (regular  wave)  solution  at  a  specified  frequency 
of  encounter,  the  overbar  denotes  the  time  average  of  one 
period,  and  ii  is  the  unit  normal  to  the  surface.  In  this  work 
we  neglect  the  effects  of  the  scattering  potential  ^g,  and  we 
consider  only  the  incident  waves  potential  This  is  a  good 

approximation  for  mostly  head  or  following  seas.  For  beam 
seas,  the  potential  ^g  is  expected  to  play  an  increasingly 
important  role.  It  should  be  mentioned,  however,  that 
periscope  depth  operations  in  beam  seas  are  not  very  common 
due  excessive  first-order  roll  motions. 

D.  NUMERICAL  SOLUTION 

Once  the  velocity  potential  for  a  distribution  of  sources 
is  formulated,  the  source  strengths  can  be  solved  for  by 
satisfying  the  body  boundary  condition.  There  are  various 
analytical  techniques  for  doing  this.  This  work  utilizes 
the  source  panel  method  developed  by  Hess  and  Smith  (Hess  & 
Smith,  1964) .  Their  method  is  detailed  below. 
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The  velocity  potential  was  derived  euad  given  in  (39)  and 


repeated  here 


<t>{P)  -^jj^G(P,Q)o{Q)dS{Q)  ,  (69) 


where  P{x,y,z)  is  the  field  or  observation  point  and 
is  the  distributed  source  point  on  the  body  surface.  We  must 
determine  a(Q)  to  satisfy  the  body  boundary  condition  on  Sjj. 
The  surface  is  represented  as  a  finite  number  of  elements  each 
having  a(Q)  constant  so  that  the  source  strength  can  be 
brought  outside  of  the  integral.  The  integral  then  depends 
only  on  the  geometry  of  the  various  panels.  With  a  finite 
number  of  unknowns  we  can  only  satisfy  the  body  boundary 
condition  at  an  equal  number  of  discrete  points. 

To  form  the  panels  the  body  surface  is  divided  into  N 
plane  quadrilaterals  with  a  constant  source  strength  over 
each  quadrilateral  for  j  «  1,2,...N.  The  (Tj  can  then  be  moved 
outside  the  integral  in  (68)  yielding  the  following  formula 
for  the  perturbation  potential 

4>IP)  G{P,Oj)dSj  ,  (70) 

where  Qj  ■  is  the  source  point  at  the  panel  and 

dSj  is  the  surface  area  of  the  panel.  The  integral  can 
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then  be  obtained  algebraically  in  terms  of  x,y,z  and  the  four 
comers  for  1  *  1,2, 3, 4.  If  (x,y,z)  is  far  from  the 

pcmel,  siirplification  can  be  used  since  the  distance  between 
the  source  point  and  the  observation  point  is  essentially 
constant  over  the  whole  panel.  The  body  boundary  conditions 
are  applied  at  one  point  in  each  quadrilateral  to  form  N 
linear  equations  in  N  unJcnowns  which  can  then  be  solved  for 
(Tj.  Once  the  source  strengths  Cj  are  determined,  the 
velocities  at  the  control  points  are  computed  from  the  formula 

tf 

Ui  =  -In  +  ^  ,  (71) 

where  VG^j  is  the  gradient  of  the  total  Green  function.  This 
includes  the  effect  of  all  three  terms  of  (36)  as  well  as  both 
sides  of  the  body.  The  Bernoulli  Equation  is  then  used  to 
give  the  dynamic  pressure 

Pi  ^  ^p(U^  -  Vi^)  .  (72) 

The  hydrostatic  pressure  is  ignored. 

In  practice  the  Hess  and  Smith  program  works  in  a  two-step 
process.  The  first  step  is  the  quadrilateral  generation. 
Points  are  input  which  specify  the  body  surface.  The  program 
then  forms  quadrilaterals  as  flat  surfaces  auid  determines  the 
normal  to  that  surface.  It  then  finds  the  null  point,  this 
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being  the  point  where  the  source  of  the  panel  produces  no 
tangential  velocity.  Finally,  it  provides  an  output  to  check 
the  points.  The  second  step  yields  the  solution.  The  program 
satisfies  the  body  boundary  condition  at  the  null  point  of 
each  quadrilateral  to  form  N  equations.  It  then  solves  for 
velocity  and  pressure  at  each  null  point  for  eui  inputted 
inflow  in  x,  y,  and  z.  Additionally,  it  will  also  obtain 
velocity  and  pressure  at  specified  off  the  body.  With  the 
pressure  at  each  panel  known,  the  drag,  lift  and  bow-up 
moments  can  be  calculated  by  integrating  over  the  appropriate 
areas  and  distances. 
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Ill 


RESX7LTS 


A.  IHTRODUCTZON 

To  inclement  the  Hess  and  Smith  source  panel  method,  a 
FORTRAN  progreua  written  by  Doctors  and  Beck  provided  the 
solution  for  all  problem  runs.  This  program  computes  the 
hydrodynamic  forces  and  moments  for  a  body  travelling  at  a 
constant  speed  by  the  Neumann- Kelvin  method.  The  key  Inputs 
to  the  program,  stored  In  a  file  IN. DAT,  Included  the 
following:  water  density,  acceleration  due  to  gravity,  body 
length,  body  beeun,  body  depth,  body  speed,  and  the  number  of 
points  longitudinally,  vertically  and  In  the  0- Integration. 
Parameters  could  be  varied  for  particular  tuns .  Additionally, 
the  body  surface  was  Input  with  a  data  file  of  longitudinal 
points,  with  a  corresponding  radius  at  each  point.  A 
subroutine  then  calculated  the  y  and  z  points  from  the  radius 
(R)  using  the  following  relations: 

y  »  R  sin(ir(IZ  -  1.0)/(NZ  -  1.0) )  (73) 


aind 

2  «  R  cos(ir(IZ  -  1.0)/(NZ  -  1.0))  -  H  ,  (74) 

where  NZ  Is  the  total  ninnber  of  vertical  points,  IZ  Is  the 
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successive  vertical  point  and  H  is  the  depth.  This  scheme 
forms  half  of  the  body  surface.  The  other  half  is  formed  by 
reflecting  the  points  about  the  vertical  centerline  in  a 
s\ibroutine.  The  program  then  utilizes  these  points  to 
generate  the  quadrilateral  pemels.  Thus,  the  Doctors  and  Beck 
program  provides  the  numerical  solution  to  the  Neumann -Kelvin 
problem  for  numerous  inputs. 

The  output  results,  stored  in  a  file  OUT. DAT,  contained 
appropriate  body  pareuneters  as  well  as  the  three  primary 
coefficients.  The  program  calculated  the  drag  and  lift  forces 
and  bow- up  moments  by  pressure  summation  over  the  body  surface 
panels.  The  moments  are  referenced  to  the  body  midpoint.  For 
the  wave  making  problem,  the  drag  coefficient,  Cp,  and  the 
lift  coefficient,  C^,  were  made  dimensionless  by  l/2pU^S, 
where  S  -  wetted  area,p  «  water  density  and  U  *  body  speed. 
The  bow-up  moment  coefficient, Cjj,  was  made  dimensionless  by 
l/2pU^SL,  where  L  ■  body  length.  For  the  incident  waves 
problem,  the  forces  and  moment  were  nondimensionalized  using 
l/2pA^(«)^S  and  l/2pA^ci)^SL,  respectively.  Here  A  ■  wave 
amplitude  and  w  «  wave  frequency.  These  nondimens ional 
coefficients  served  as  the  basis  for  evaluating  the  effects  on 
the  body. 

The  model  used  to  generate  the  input  data  points  was  a 
submarine  model,  SUBOFF.  Developed  by  the  Defense  Advanced 
Research  Projects  Agency  (DARPA) ,  the  model  was  originally 
developed  to  evaluate  various  flow  field  predictions  for  an 
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axisymmetric  hull,  both  with  and  without  appendages.  It  was 
intended  to  con^are  predictions  with  model  experimental  data. 
For  all  runs  in  this  work,  data  for  a  bare  hull  (no 
appendages)  model  were  utilized.  The  key  parameters  of  this 
model  are  a  length  of  14.2917  ft,  a  1.6667  ft  midbody  diameter 
and  a  displacement  of  24.692  ft^.  To  model  the  hull  surface, 
various  station  points  (longitudinal  points)  were  provided 
along  with  their  corresponding  radius  (fraction  of  midbody 
diameter)  .  These  points  are  tedjulated  in  Table  I.  These 
station  points  were  converted  into  corresponding  x  values, 
with  station  0.0000  at  the  bow  (+L/2)  and  station  20.4167  at 
the  stem  (-L/2)  .  The  radius  values  were  used  to  produce  the 
y  and  z  values.  A  drawing  of  the  DARPA  SUBOFF  model  is 
included  as  Figure  2.  The  remaining  sections  of  this  chapter 
detail  the  convergence,  the  wave  making  problem  and  the 
incident  waves  problem  results.  The  convergence  section 
esqplores  the  optimum  number  of  vertical  and  6  steps  required 
to  converge  to  an  acceptable  solution.  The  wave  making 
problem  section  examines  the  depth  dependence  of  each 
coefficient  over  a  range  of  speeds.  Finally,  the  incident 
waves  section  looks  at  two  areas.  The  first  is  the  effect  on 
the  coefficients  when  the  wave  period  is  varied.  The  second 
area  shows  the  effect  on  the  coefficients  when  the  wave 
frequency  of  encounter  and  the  wave  direction  is  altered. 


TABLE  I.  SUBOPF  STATION  POINTS  AND  RADII 


STATION  POINTS 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


1.0 


2.0 


3.0 


4.0 


7.7143 


10.0 


15.1429 


16.0 


17.0 


18.0 


19.0 


20.0 


20.1 


20.2 


20.3 


20.4 


-2(L..41£2. 


RADIUS  FRACTION 


0.00000 


0.29058 


0.39396 


0.46600 


0.52147 


0.56627 


0.60352 


0.63514 


0.70744 


0.84713 


0.94066 


0.99282 


1.00000 


1.00000 


1.00000 


0.97598 


0.81910 


0.55025 


0.26835 


0.11724 


0.11243 


0.10074 


0.07920 


0.03178 


O.OOOOQ 
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(u)  i»tP"a 


I 


Figure  2:  Drawing  of  the  DARPA  SUBOFF 
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B.  CONVERGENCE 


The  first  step  in  producing  accurate  results  requires  a 
consistent  number  of  integration  increments  in  the  vertical 
and  directions  be  used  for  each  run.  The  longitudinal  or  x 
increments  (NX)  are  fixed  by  the  station  points  of  the  DARPA 
SUBOFF  model  at  25.  A  series  of  coir^uter  runs,  using  steuidard 
input  data,  was  conducted  to  determine  the  number  of  vertical 
steps,  NZ,  euid  the  nxunber  of  9 -steps,  N6.  The  input 
conditions  were  standardized  as  follows: 

•  p  ■  1.94  slugs/ft^ 

•  g  -  32.2  ft/sec^ 

•  L  -  14.3  ft 

•  D  -  1.67  ft 

•  H  -  0.15L  -  2.15  ft 

•  U  -  10.0  ft/sec 

•  NX  -  25 

The  convergence  for  the  steps  (N0)  was  first  conducted. 
The  depth- to- submergence  ratio,  H/L,  of  0.15  was  chosen  to 
closely  correspond  to  the  value  of  0.16  that  Doctors  and  Beck 
utilized  in  their  study  of  a  submerged  spheroid.  The  spheroid 
had  a  depth- to- length  ratio  of  0.2,  whereas  the  SUBOFF  has  a 
ratio  of  0.12.  Five  runs  were  conducted,  using  a  fixed  value 
of  vertical  steps,  NZ  -  8.  The  runs  varied  NH  from  4  to  20, 
in  increments  of  4.  Figures  3  to  5  clearly  show  a  large 
change  in  all  three  coefficients  for  values  of  NH  below  8. 
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Cougar ing  the  percent  change  between  NB  *  4  emd  NB  -  8  for 
each  plot,  the  drag  coefficient  (C^)  varied  by  only  0.80%,  the 
lift  coefficient  (Cj,)  changed  by  11.9%  and  the  moment 
coefficient  {Cj^)  rose  by  1.44%.  Above  NB  -  8,  the  plots  show 
a  nearly  constant  value.  steadied  at  -2.76E-2,  leveled 
at  1.51E-3  and  Cq  stabilized  at  7.84E>3.  Since  a  value  of  NB 
«  16  yielded  effectively  convergent  results  for  each 
coefficient  (to  three  significant  figures) ,  this  value  was 
chosen  as  the  standard  for  all  subsequent  runs.  This 
selection  was  important,  as  the  computational  time  varied  with 
the  square  of  NB. 

To  study  the  effects  of  the  vertical  steps,  NZ,  the  value 
of  NB  was  fixed  at  16.  NZ  was  varied  from  4  to  9  in  steps  of 
one.  Figures  6  through  8  show  the  convergence  plots.  In  each 
case,  the  plot  has  not  leveled  off  at  NZ  «  9.  However,  the 
slope  in  each  plot  has  significantly  decreased  by  NZ  «  8  from 
its  initial  value  at  NZ  «  4.  Drag  slope  dropped  from  8.47E>4 
to  1.14E-4,  lift  slope  changed  from  2.25E-4  to  3.03E-5,  and 
moment  slope  varied  from  -2.99E-4  to  -3. HE- 5.  Between  NZ  « 
8  and  NZ  «  9,  the  values  of  the  drag,  lift  and  moment 
coefficients  vary  by  only  1.45%,  2.00%,  and  1.13%, 
respectively.  A  value  of  NZ  >  8  was  chosen  as  the  standard 
for  future  runs  based  on  this  reasonable  convergence. 
Choosing  more  vertical  points  would  not  significantly  improve 
the  solution.  Additionally,  it  would  also  increase  the 
con^utational  time  since  it  is  also  proportional  to  the  square 
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of  the  number  of  panels.  Also,  the  ratio  of  NX/NZ  -  25/8  ■ 
3.125  is  smaller  than  the  ratio  of  4  used  in  the  Doctors  and 
Beck  study  to  achieve  reasoncible  convergence.  The  smaller 
ratio  is  better  since  the  quadrilaterals  that  are  formed  will 
have  a  more  uniform  shape  than  a  body  with  a  higher  aspect 
ratio. 

Values  of  N0  •  16  and  NZ  -  8  were  selected  as  the 
standards  for  subsequent  runs.  As  this  section  detailed, 
these  values  give  good  convergence  results  in  reasonable 
computational  time.  Using  these  values,  the  runs  to  determine 
the  wave  making  effects  of  the  body,  were  performed. 
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DRAG  COGI-|-|Cll£NT 


Figure  3:  Theta  Convergence  for  Drag  Coefficient 
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llf-T  COEFFICIfENT 


TriETA  CONVERGENCE  FOR  UFT 


4  6  3  10  12  U  16  18  20 

INTEGRATiCN  STEPS  IN  THETA  (NH) 

Figure  4:  Theta  Convergence  for  Lift  Coefficient 


MOMENT  COEFFICIENT 


DRAG  COlEri-lClENT 


2  CONVEr.GSNCc  FOR  DRAG 


Figure  6:  Z  Convergence  for  Drag  Coefficient 
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X  10’^ 


Z  CONVERGENCE  FOR  LIFT 


MOMENT  COEFFICIENT 


X  10' 


2  CONVERGENCE  FOR  MCMENI 


Figure  8:  Z  Convergence  for  Moment  Coefficient 
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C.  WAVE  MA1CZW6  PROBLEM 


Having  established  N6  and  NZ,  the  DARPA  SUBOFF  model  is 
input  into  the  program  to  solve  the  wave  making  problem.  The 
drag,  lift  and  moment  coefficients  are  calculated  as  a 
function  of  the  body  speed  (Froude  number)  for  varying  depths. 
The  Froude  number  is 

F  - 

J  • 


where  U  is  the  body  speed  in  ft/sec,  g  is  the  acceleration  due 
to  gravity  and  L  is  the  body  length.  The  depths  (H)  were 
chosen  to  be  a  fraction  of  the  ship  length,  ranging  from  one 
tenth  to  one  half.  The  following  depth  ratios  (H/L)  were 
chosen  to  give  a  good  spread  of  values,  with  more  iruns 
performed  at  shallower  depths: 

•  0.100 

•  0.110 

•  0.120 

•  0.125 

•  0.135 

•  0.150 

•  0.175 

•  0.200 

•  0.225 
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•  0.250 


•  0.300 

•  0.400 

•  0.500 


Figures  9,  10,  and  11,  plot  the  drag,  lift  and  moment 
coefficients  versus  Froude  number.  The  shallowest  (0.1)  and 
the  deepest  (0.5)  depth  ratios  are  ladjelled,  with  each 
successive  curve  being  the  sequential  plot  from  the  above 
list.  The  Froude  number  is  plotted  from  0.18  to  0.75  which 
corresponds  to  a  range  of  4  to  16  ft/sec.  A  spline  technique 
was  utilized  to  smooth  the  curve. 

The  first  plot.  Figure  9,  graphs  the  Drag  Coefficient 
versus  Froude  Number.  The  drag  coefficient  indicates  the 
force  which  is  opposing  the  body  in  the  longitudinal 
direction.  The  oscillatory  nature  of  the  curves  is  readily 
apparent  especially  at  the  shallower  depths.  At  deeper 
depths,  the  near  surface  effects  are  greatly  diminished  and  it 
exhibits  small  oscillations  about  a  steady  value  of  edsout 
0.005.  Two  pealcs  are  clearly  seen  at  about  F  -  0.3  and  F  - 
0.5.  The  magnitude  of  the  peaks  decreases  and  the  peak  shifts 
to  the  right  as  depth  is  increased.  With  the  exception  of  the 
shallowest  depth,  the  second  peak  (higher  speeds)  is  of  a 
greater  magnitude  than  the  first  pecdc.  For  the  shallowest 
depth,  the  ratio  of  the  highest  to  lowest  coefficient  is  about 
3.  This  ratio  drops  off  to  just  over  1  for  the  deepest  depth. 
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Also  noteworthy  is  the  observation  that  at  F  «  0.2  and  P  » 
0.35,  a  minimum  coefficient  is  reached,  regardless  of  depth. 
Clearly,  the  drag  coefficient  markedly  decreases  with  depth, 
since  the  Interaction  with  the  calm  free  surface  exponentially 
decays . 

A  similar  oscillatory  behavior  is  observed  in  Figure  10, 
Lift  Coefficient  versus  Froude  Number.  At  shallow  depths  and 
low  speeds,  a  large  lift  is  generated.  The  positive  values  of 
lift  represent  an  upwards,  "sea- suction"  force  on  the  body. 
It  is  this  force  which  actual  submarines  must  counter  to 
remain  submerged  at  periscope  depth.  Two  local  maxima  and  two 
local  minima  are  observed.  As  with  the  drag  coefficient,  the 
maocima  shift  to  the  higher  speeds  at  deeper  depths.  However, 
in  contrast  with  the  drag  coefficient,  the  lift  significantly 
decreases  with  increasing  speed  and  becomes  negative  (downward 
force)  at  speeds  aUoove  F  »  0.55.  For  the  shallowest  depth, 
the  second  local  maximum  is  almost  four  times  smaller  than  the 
first  local  maximum.  This  correlates  well  to  actual 
siibmarines ,  in  which  it  is  easier  to  remain  submerged  at 
higher  speeds.  Again,  at  the  deepest  depth,  the  oscillations 
appear  as  a  straight  line  using  the  scale  in  Figure  10.  Its 
lowest  coefficient  is  0.0462E-3  at  the  lowest  speed,  and  it 
reaches  a  maximum  of  0.126E-3  at  F  -  0.466  (a  ratio  of  2.74) . 
So,  the  deepest  depths  do  oscillate  around  a  fairly  constant 
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value,  but  the  effect  on  the  body  Is  negligible  compared  with 
shallower  depths.  Thus,  the  shallowest  depths  and  the  lowest 
speeds  exert  the  greatest  lift  force  on  the  body. 

The  Moment  Coefficient  versus  Froude  Number,  Figure  11, 
likewise  displays  an  oscillatory  characteristic.  The  negative 
values  for  the  coefficient  Indicate  that  a  bow- down  moment  is 
Induced  on  the  body,  rotating  the  stem  end  towards  the 
surface.  This,  too,  agrees  with  the  motion  of  actual 
submarines  first  broaching  the  surface  at  the  stem.  Two 
local  maxima  and  two  local  minima  are  present  as  with  the 
other  two  coefficients.  Also,  the  shallowest  depths  produce 
the  largest  fluctuation  in  the  coefficient.  However,  the 
magnitude  of  the  coefficient  greatly  increases  as  the  speed  is 
increased.  The  oscillations  appear  less  distinct  as  the  depth 
Increases,  but  the  magnitude  is  clearly  greater  for  higher 
speeds.  The  ratio  of  the  local  maxima  for  the  shallowest 
depths  is  about  three  to  one  and  drops  to  just  over  one  for 
the  deepest  depths.  As  seen  with  the  other  two  coefficients, 
the  local  maxima  shift  to  the  right  (higher  speeds)  with 
increasing  depth.  In  other  words,  at  shallower  depths  the 
larger  moment  coefficient  is  produced  at  lower  speeds.  Also 
of  note  is  the  greater  moment  coefficient  magnitude  generated 
at  lower  speeds  (below  F  -  0.35)  as  depth  is  increased  2d>ove 
0.15L.  The  moment  coefficient  causes  the  body  to  pitch  down 
for  all  speeds  and  depths. 
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The  wave  making  problem  plots  clearly  display  the  results 
of  summing  the  pressure  distributions  around  the  body.  All 
three  coefficients  are  oscillatory  in  nature.  As  the  body  is 
placed  in  deeper  depths,  the  coefficients  become  smaller  in 
magnitude.  As  speed  is  increased,  drag  and  the  bow-down 
moment  increase,  while  lift  decreases.  The  results  correspond 
well  to  the  operating  experience  of  actual  submarines  at 
shallow  depths. 
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Wave  Making:  Drag  Coefficient  versus  Froude 
Number 
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^  ^q-sWAVE  MAKING:  MOMENT  vs  FECUDE  NUMBER  (VARIABLE  DEPTH) 


S.  INCIDENT  NAVES  PROBLEM 


After  examining  the  body  vrave  mciking  problem,  the  effect 
of  incident  waves  on  the  body  is  explored.  The  SUBOFF  model 
was  input  to  provide  the  program  the  body  shape  with  which  the 
incident  waves  interacted.  Keeping  only  the  second  term  of 
equation  (65)  and  sxibstituting  the  encounter  frequency  from 
equation  (67)  yielded  the  following  pressure  expression 

p  =  .  (76) 


This  pressure  distribution  was  summed  over  all  the  panels  to 
produce  the  drag,  lift  and  moment  coefficients.  The  drag  and 
lift  coefficient  were  normalized  by  l/2p(i>^A^S  and  the  moment 
coefficient  was  normalized  by  l/2p(i)^A^SL.  For  the  various 
runs  in  this  section,  the  depth,  H,  the  body  speed,  U,  the 
wave  amplitude.  A,  the  wave  direction,  6,  and  the 
nondimens ional  wavelength,  X/L,  were  input  to  the  program. 
Aitplitude,  A,  was  fixed  at  1.0  and  X/L  varied  from  0.5  to  2.5 
for  all  runs.  The  wavelength  was  directly  obtained  from  the 
nondimens  ional  wavelength  and  the  wavenumber,  k  «  2ir/X, 
followed.  The  deep  water  wave  frequency  »  k/g  was  then 
obtained  and  used  to  calculate  the  encounter  frequency, 

-  (  «^/g)UcosC,  where  0  -  0®  is  following  seas  and  0  - 
180®  is  head  seas.  These  parameters  determined  the  pressure 
distribution  through  (75) .  All  of  the  plots  in  this  section 
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graph  the  coefficients  versus  the  nondimens ional  wavelength, 
X/L.  Speeds  are  expressed  as  a  Froude  number,  depth  as  a 
fraction  of  overall  length  euid  the  wave  direction,  6,  as  an 
angle  in  degrees.  A  positive  lift  coefficient  represents  an 
upwards  (sea  suction)  force  and  a  positive  moment  coefficient 
shows  a  bow-up  moment.  Three  sets  present  the  effects  on  the 
three  coefficients.  The  first  set  examines  the  effect  of 
depth  as  speed  and  wave  direction  are  held  constant.  The  next 
set  varies  the  wave  direction  at  a  constant  speed  and  depth, 
while  the  last  set  changes  the  body  speed  with  depth  and  wave 
direction  being  fixed. 

The  first  incident  wave  set  looks  at  the  exponential 
nature  of  the  depth  term  for  a  given  set  of  conditions.  Hull 
speed  was  fixed  at  F  «  0.2,  and  $  set  at  90**  (beam  seas). 
This  removed  the  effect  of  the  body  speed  from  the  frequency 
of  encounter.  Depth  was  then  varied  between  a  tenth  and  a 
half  of  the  overall  body  length.  Figures  12  through  14 
present  the  Drag,  Lift  and  Moment  Coefficients  versus 
Nondimensional  Wavelength,  respectively.  In  each  plot  the 
coefficients  rapidly  decay  as  the  depth  is  increased.  The 
coefficients  are  negligible  at  a  depth  greater  than  0.5X  (0.5L 
times  l.OX/L).  This  validates  the  statement  made  in  the 
problem  formulation  section  that  there  is  a  negligible  motion 
beneath  a  depth  of  0.5X.  As  the  wavelength  increases,  the 
coefficients  also  increase  at  the  deepest  depths,  showing  that 
the  longer  waves  have  a  larger  effect  on  the  body.  However, 
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at  the  shallowest  depth  the  coefficients  decrease  with 
increasing  wavelength.  Both  the  lift  cuid  moment  coefficients 
peak  at  about  1.2X/L  then  steadily  drop  as  the  waves  lengthen. 
Overall,  these  plots  clearly  show  the  exponential  decay  of  the 
coefficients  with  depth. 

Figures  15  to  17  show  the  effects  on  the  coefficients  when 
the  wave  direction  is  varied.  For  these  runs,  depth  was  fixed 
at  one  tenth  of  overall  length  and  body  speed  was  set  at  F  « 
0.2.  Wave  direction  was  varied  from  0”  (following  seas)  to 
180**  (head  seas)  in  45^  increments  and  are  labelled  as  such. 
In  each  plot,  the  9  «  090  corresponds  to  the  H  -  O.IL  curves 
of  the  varioLble  depth  plots.  In  the  Drag  Coefficient  versus 
Nondimens ional  Wavelength  plot  (Figure  15) ,  maximum  drag  is 
produced  with  head  seas  and  steadily  drops  off  as  the  seas 
move  aft.  The  following  seas  are  an  order  of  magnitude 
smaller  than  the  head  seas.  As  the  wavelength  increases,  the 
drag  drops  as  the  waves  appear  longer  to  the  body.  Figures  16 
cuid  17,  the  Lift  and  Moment  Coefficients  versus  Nondimens  ional 
Wavelength,  show  identical  trends  to  one  another.  For  head 
seas,  both  plots  peak  around  a  wavelength  of  one  then  steadily 
decrease  as  the  wavelength  increases.  The  head  seas  produce 
over  twice  the  magnitude  of  the  beam  seas.  This  correlates 
well  to  at  sea  experience  of  being  able  to  better  maintain 
depth  with  beam  seas  than  with  head  seas.  For  the  following 
seas,  the  lift  and  moment  coefficients  both  slowly  Increase 
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with  wavelength,  but  lie  significantly  below  the  beam  seas 
curve.  Thus,  head  seas  produce  the  largest  coefficients  and 
following  seas  yield  the  smallest  coefficients. 

Next,  the  effects  on  varying  the  body  speed  were  examined. 
Figures  18  to  20  plot  the  Drag,  Lift  and  Moment  Coefficients 
versus  Nondimens ional  Wavelength.  In  these  runs,  depth  was 
fixed  at  one  tenth  of  overall  length  and  9  set  to  180<*  (head 
seas).  Speed  was  varied  from  F  -  0.0  to  0.8  in  0.2 
increments.  The  F  «  0.2  curves  correspond  to  the  9  *  180 
curves  on  the  variable  wave  direction  plots  (Figures  15  to 
17).  Also,  the  F  -  0.0  curves  correspond  to  the  H  «  O.IL 
curves  on  the  variable  depth  plots  (Figures  12  to  14) ,  since 
the  effect  of  hull  speed  was  taken  out  of  these  by  using  becim 
seas.  Each  plot  clearly  shows  that  the  greater  speeds  produce 
the  largest  coefficients.  From  the  lowest  to  the  highest 
speeds,  the  coefficients  change  by  an  order  of  magnitude. 
Additionally,  as  the  wavelength  is  raised,  each  of  the 
coefficients  slowly  decays  in  magnitude.  Thus,  in  examining 
the  incident  waves  effects  alone,  a  higher  hull  speed  or  a 
shorter  wavelength  would  cause  the  most  drag,  the  greatest 
lift  and  the  largest  bow-down  moment. 

The  results  of  the  incident  wave  runs  clearly  showed  the 
effects  on  the  SUBOFF  model  for  various  input  conditions.  The 
coefficients  all  dropped  exponentially  with  depth  as  was 
predicted.  Below  0.5X,  the  incident  waves  are  negligible. 
The  head  seas  produced  the  greatest  coefficients  at  all 
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wavelengths.  Seas  aft  of  the  beam  are  significantly  less 
important.  The  higher  hull  speeds  yielded  the  largest 
coefficients  at  each  wavelength.  The  lowest  speed  runs  were 
equivalent  to  the  beam  seas  plots. 
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ORAG  COEFI-ICIEN  r 


NONDIMENSICNAL  LAMBDA 


Figure  12:  Incident  Wave:  Drag  Coefficient  versus 

Nondimens ional  Wavelength  (Variable  Depth) 
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INC  WAVE;  UFTvs  LAMBDA  O/AP.IABLE  DEPTH) 


Figure  13:  Incident  Wave:  Lift  Coefficient  versus 

Nondimensional  Wavelength  (Variable  Depth) 


MOME 


INC  WAVE;  MOMENT  vs  LAMBDA  (VARIABLE  DEPTH) 


NONDIMENSiCNAL  LAMBDA 


Figure  14:  Incident  Wave:  Moment  Coefficient  versus 

Nondimens ional  Wavelength  (Varicdjle  Depth) 


INC  WAVE;  DRAG  vs  LAMBDA  0/AP.IABLS  WAVE  DIRECTION) 


I IFT  COEFFICIENT 


INC  WAVE;  UFTvs  LAMBDA  0/ARIABLH  WAVE  DIRECTION) 


NCNDIMENSICNAL  LAMBDA 


Figure  16:  Incident  Wave:  Lift  Coefficient  versus 

Nondimens ional  Wavelength  (Variable  Wave 
Direction) 
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^  ^q-3  INC  WAVE:  MOMENT  vs  LAMBDA  (VARIABLE  WAVE  DIRECTION) 


0.5  1  1.3  2  2.5 

NONDIMENSIONAL  L^^iMBDA 


Figure  17:  Incident  Wave;  Moment  Coefficient  versus 
Nondimens ional  Wavelength  (Variable  Wave 
Direction) 
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INC  WAVE;  CRAG  vs  LAMBDA  (VARIABLE  SPEED) 


Figure  18:  Incident  Wave:  Drag  Coefficient  versus 

Nondimensional  Wavelength  (VaricQjle  Speed) 
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I.IFT  COEFFICIENT 


INC  WAVE:  UFT  vs  LAMBDA  (VARIABLE  SPEED) 


Figure  19:  Incident  Wave:  Lift  Coefficient  versus 

Nondimens ional  Wavelength  (Variable  Speed) 
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MOMENT  COEFFICIENT 


INC  WAVE:  MOMENT  vs  LAMBDA  (VARIABLE  SPEED) 


Figure  20:  Incident  Wave:  Moment  Coefficient  versus 

Nondimens ional  Wavelength  (Variaible  Speed) 
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E.  smouotY 


The  results  chapter  constitutes  the  major  part  of  this 
work.  Using  the  DARPA  SUBOFF  model  of  a  bare  submarine  hull, 
the  hydrodynamic  forces  and  moments  produced  from  various 
sources  was  examined.  A  set  of  convergence  runs  established 
the  number  of  vertical  and  0- integration  points  at  8  and  16, 
respectively.  These  were  utilized  for  the  subsequent  wave¬ 
making  runs  and  the  incident  wave  runs. 

The  wave -making  runs  produced  oscillatory  plots  for  each 
coefficient.  The  drag  coefficient  represented  the 
longitudinal  force  opposing  the  hull's  motion  and  this 
steadily  declined  with  increasing  depth.  The  lift  coefficient 
showed  the  upwards,  sea- suction  force  present  on  the  body.  At 
the  shallowest  depths  and  the  slowest  speeds,  these  sea- 
suction  forces  were  the  greatest.  The  moment  coefficient 
explained  how  much  the  body  would  pitch  up  or  down.  As  speed 
increased,  the  bow- down  moment  increased.  Of  these  effects, 
the  lift  force  has  the  greatest  bearing  on  a  submarine 
maintaining  depth.  Thus,  using  the  greatest  speed  is  the  most 
advantageous  to  prevent  broaching  the  submarine. 

Finally,  the  incident  wave  znns  showed  the  effects  on  the 
body  from  incident  waves  on  an  otherwise  calm  free  surface. 
The  incident  waves  are  e3q>onentially  attenuated  as  depth  is 
increased  and  at  depths  greater  than  half  the  wavelength,  the 
effects  are  negligible.  Head  seas  produce  a  larger  lift  and 


69 


moment  them  do  beeun  or  following  seas,  maJelng  it  more 
difficult  to  maintain  depth  when  driving  a  submarine  into  the 
seas.  Upon  excunlning  the  role  of  hull  speed,  it  was  found 
that  the  greatest  speeds  produced  the  most  drag,  the  highest 
lift  and  the  largest  moment.  Thus,  the  higher  hull  speed 
would  increase  the  severity  of  the  incident  waves'  lift. 

Combining  the  wave  making  and  incident  wave  results  yields 
the  total  response  on  the  hull.  The  total  drag  on  the  body 
reflects  the  resistance  to  the  body's  forward  motion.  The 
wave  making  runs  produced  positive  (opposing)  forces  on  the 
body  while  the  incident  wave  runs  yielded  very  small  negative 
(aiding)  forces  on  the  body.  The  drags  for  each  were 
calculated  to  make  realistic  comparisons.  For  the  wave  making 
runs  the  drag  (D^^)  is 

-  Cj^^pU^S  ,  (77) 

and  for  the  incident  waves  the  drag  (Dji,)  is 

*  Cijlpu^A^S  .  (78) 

For  the  low  speeds  (F  -  0.2)  with  head  seas,  -  0.765  Ibf 
cmd  Djf,  «  -0.00928  Ibf.  The  two  forces  are  of  opposite  signs, 
but  the  magnitude  of  the  wave  making  drag  greatly  outweighs 
the  incident  wave  drag.  A  similar  effect  is  also  seen  at 
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higher  speeds  (F  *  0.7)  with  -  1.61  Ibf  and  Dj,,  -  -0.0325 
Ibf .  Thus,  the  majority  of  the  drag  on  the  submerged  body  is 
from  wave  making  effects. 

The  total  moment  produced  on  the  body  directly  affects  the 
ability  of  the  siibmarine  to  keep  the  stem  (mainly)  or  bow 
from  breaking  the  free  surface.  The  wave  making  mns  produced 
negative  (bow- down)  moments  and  the  incident  waves  yielded 
positive  (bow-up)  moments.  The  moments  for  each  run  were 
calculated  to  make  relevant  comparisons.  For  the  wave  making 
runs  the  moment  (M^^)  is 

~  •  (79) 

cuid  for  the  incident  waves  the  moment  (Mj^,)  is 

Wjw  =  .  (80) 

For  the  low  speeds  (F  -  0.2)  with  head  seas,  ■  -1.44  ft- 
Ibf  and  >  2.10  ft- Ibf.  This  shows  that  the  two  moments 
are  roughly  equivalent  and  of  opposite  sign,  thus  cancelling 
out  one  another.  As  the  seas  move  aft,  drops  to  0.93  ft- 
Ibf  off  the  beam  and  falls  to  0.22  ft-lbf  for  following  seas. 
For  these  cases,  a  total  bow-down  moment  is  produced,  forcing 
the  stem  up  to  the  surface.  As  speed  is  raised  to  F  -  0.7, 
all  the  moments  increase,  except  for  the  beam  incident  waves 
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which  are  Independent  of  speed.  rises  dramatically  to  - 
123.4  £t-lb£  while  only  increases  to  7.30  £t-lb£  for  head 
seas  and  to  0.53  ft-lbf  for  following  seas.  Thus,  a  much 
greater  bow- down  moment  is  produced  at  higher  speeds, 
resulting  in  more  submarine  control  surfaces  and  varicible 
ballast  being  needed  to  prevent  the  stem  from  breaking  the 
free  surface. 

Even  of  greater  in^ortance  to  submarine  depth  control  is 
the  effect  of  the  combined  lift  on  the  body.  This  is  the 
primary  item  of  interest  since  this  upward  force  must  be 
overcome  to  prevent  the  entire  body  from  broaching  the  free 
surface.  A  combination  of  ship's  speed,  control  surfaces  and 
variable  ballast  are  utilized  to  counter  the  lift  force.  In 
the  wave  making  runs,  higher  ship  speeds  produced  less  lift, 
while  greater  speeds  yielded  more  lift  in  the  incident  wave 
runs.  To  make  the  results  meaningful,  the  lift  forces  in  each 
case  must  be  compared.  For  the  wave  making  runs,  the  lift 

,  (81) 


and  for  the  incident  wave  runs  the  lift  is 

^ jw  “  C'x, p  A  ^  S 


(82) 


For  the  case  of  head  seas  with  body  speed  being  F  «  0.2  and 
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nondimensional  wavelength  of  one,  hj^  «  27.64  Ibf  and  - 
3.02  Ibf.  When  the  seas  come  from  the  beam,  Ljf,  -  12.37  Ibf 
euid  when  the  seas  are  following,  Lj|f  «  3.09  Ibf.  In  all 
cases,  both  lifts  are  positive  eind  the  body  is  drawn  to  the 
surface.  When  speed  is  increased  to  F  •  0.7  in  head  seas, 
significantly  increases  to  92.80  Ibf,  but  L|q(  changes  sign  cuid 
is  of  the  same  order  of  magnitude,  -4.93  Ibf.  The  upwards 
force  from  the  incident  head  waves  will  greatly  overcome  the 
negative  lift  generated  by  the  body.  But,  as  the  sea 
direction  moves  aft,  the  lift  generated  by  the  incident  waves 
drastically  drops.  For  beam  seas,  Lji,  remains  constant 
(independent  of  speed)  and  for  following  seas  Lj|f  rises  to 
€.87  Ibf.  Thus,  sea  direction  plays  a  major  role  in 

determining  the  lift  which  the  submarine  must  counteract. 
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TV.  PARAMETRIC  STUDIES 


A.  INTRODUCTION 

The  previous  chapter  detailed  the  results  of  the  wave 
making  problem  euid  the  incident  waves  problem  for  a  given 
body,  namely  the  DARPA  SUBOFF  model.  This  chapter  explores 
the  effects  of  changing  the  body  size  and  shape.  The  body  is 
divided  into  three  sections,  with  each  section  being  defined 
by  a  set  of  coefficients.  By  varying  the  overall  length  in 
one  case  and  by  changing  the  maximum  dleimeter  in  a  second 
case,  the  effects  on  the  drag,  lift  and  moment  coefficients 
are  examined. 

The  body  of  revolution  hull  form  was  selected  to  model  a 
submarine  since  it  is  described  by  simple  geometric  forms 
which  are  developed  from  elementary  mathematical  equations. 
Essentially,  all  cross  sections  of  the  hull  are  circular.  The 
development  of  this  approach  is  articulated  in  an  article  on 
submarine  design  concepts  (Jackson,  1992) .  The  hull  form  is 
composed  of  three  sections.  The  forward  end  called  the 
entrance,  the  parallel  middle  body  and  the  after  end  called 
the  run.  The  entrance  is  a  portion  of  a  parabola  of 
revolution  with  length,  Lf,  of  2.4  dicuneters.  The  run  is  a 
portion  of  an  ellipsoid  of  revolution  with  length,  of  3.6 
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diameters.  The  parallel  middle  body  is  a  cylinder  of  meocimum 
diameter,  O,  with  length,  Lp^jg-The  sum  of  Lf  and  Lpjjg  is 
the  overall  length,  L. 

If  one  were  to  use  equations  for  true  ellipsoids  and 
parabolas,  the  entrance  and  run  would  be  too  fine  for  a  modem 
submarine  hull  form.  The  body  size  is  increased  by  using 
large  exponents  (n^  and  n,^)  in  the  equations  below. 

<W) 

Z  Jjf 

y,  .  |[1  -  (84) 


Figure  21  is  a  drawing  of  the  hull  shape  which  details  the 
various  lengths.  Here  and  x^  are  the  offsets  from  the 
maximum  diameter,  and  y^  and  are  the  hull  radii  at  their 
respective  offsets. 


I 

I 


I 

Figure  21:  Submarine  Hull  Shape 
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Using  these  parameters,  the  total  volxime  of  the  body  can  be 
calculated  as  the  sum  of  the  volxunes  of  the  entrance,  Vf,  the 
run,  and  the  parallel  middle  body,  Vpjjg: 

Vf  ‘  Tr{^)^  [Cpf*2.A*D]  (85) 

V,  «  IT  [Cpa*3.6*D]  (86) 


VpMB  =  -  6*D] 


(87) 


where  Cpf  and  Cp^  are  the  prismatic  coefficients  for  the 
entrance  and  run,  respectively.  Sxibstituting  y^  and  y^  for 
D/2  and  integrating  over  the  length  of  each  section  yields 


Cpf  =  f(l  -  (x)”^)  ^  dx 


'pa 


jd  -  (Xi)”-)^  dxi 


(88) 


(89) 
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where  x  »  ^Cf/Lf  and  x^  «  x^/L^.  Adding  the  three  volumes  and 
rearranging  results  in 

V  «  t3.6Cpa  +  — ^  +  2.4Cpf]  .  (90) 

In  performing  the  parametric  studies,  the  volume 
(displacement)  was  kept  constant  said  set  equal  to  the 
SUBOFF's  displacement.  The  SUBOFF  model  was  slightly 
different  than  the  bodies  of  revolution  described  above,  in 
that  the  entrance  was  twice  the  diameter  and  the  run  was  2.18 
times  the  diameter.  Yet,  the  shape  of  each  SUBOFF  section  was 
certainly  circular.  In  the  first  study,  the  displacement  and 
the  diameter  of  the  body  were  fixed  and  the  length  was  varied 
through  the  use  of  different  exponents  for  the  entrance  and 
the  run.  In  the  second  study,  the  displacement  and  the 
overall  length  were  fixed  and  the  diameter  changed  with 
various  exponents.  In  both  studies,  wave  making  and  incident 
wave  runs  were  conducted. 

B.  EFFECTS  OF  LENGTH 

The  first  parametric  study  examined  the  effects  on  the 
drag,  lift  and  bow-up  moment  coefficients  as  the  overall 
length  was  varied.  The  displacement  was  fixed  at  24.692  ft^ 
and  the  maximum  diameter  held  at  1.667  ft.  All  parameters 
remained  the  same,  except  that  overall  length  varied  in  each 
run  and  depth  was  set  at  one  tenth  of  the  length.  For  each 
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run,  equal  values  of  the  exponents  n^  and  n^  were  chosen. 
From  these  exponents,  the  prismatic  coefficients  were  obtained 
from  (68)  amd  (89).  Then,  from  (90),  the  length  of  the 
parallel  middle  body,  Lpj^g,  was  found  which  set  the  overall 
length,  L.  Table  II  summarizes  the  exponents,  prismatic 
coefficients  amd  overall  lengths  utilized  in  this  parametric 
study.  As  the  exponent  increased  in  value,  the  entrance  and 
the  run  became  wider  and  rounder;  i.e.,  more  cylindrical 
throughout.  Also,  the  overall  body  length  became  shorter. 
The  longitudinal  stations  and  their  respective  radii  were 
calculated  and  converted  into  the  correct  input  format.  The 
y  and  z  values  were  calculated  as  in  previous  runs.  As 
before,  the  computer  program  output  the  three  coefficients 
which  are  shown  in  Figures  22  to  33.  The  graphs  plot  the 
drag,  lift,  and  moment  coefficients  versus  speed  (Froude 
niimber)  for  the  wave  making  runs  and  plot  the  coefficients 
versus  nondimens ional  wavelength  for  the  incident  wave  runs. 
In  each  plot,  the  n^  -  n£  -  2.0  and  the  n^  -  n£  -  4.0  curves 
are  labelled,  with  the  other  four  curves  lying  sequentially 
between  these  two  labelled  curves.  Each  plot  will  be 
separately  detailed. 
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TABLE  II.  PRISMATIC  COEFFICIENTS  AND  OVERALL  LENGTH 


1 

Cpa 

Cpf 

L 

1  2.0 

0.5333 

0.6667 

15.5  1 

1 

0.5952 

0.7500 

14.7  1 

1 

0.6429 

0.8061 

14.2 

1 

0.6806 

0.8455 

13.9 

4.0 

0.7111 

0.8740 

13.6 

4.5 

0.7366 

0.8944 

13.3 

Figures  22  to  24  detail  the  wave  making  runs.  The  Wave 
Making  Drag  Coefficient  versus  Froude  Number  plot  (Figure  22) , 
shows  an  oscillatory  character.  As  the  exponents  are  raised, 
the  curves  show  larger  fluctuations  in  amplitude. 
Additionally,  the  local  minima  and  maxima  shift  to  the  right 
as  the  exponents  increase.  The  general  shape  of  each  curve  is 
identical,  with  each  curve  slowly  decreasing  in  magnitude 
ed>ove  cQ30ut  F  >  0.5.  Overall,  the  magnitude  increases  by  a 
factor  of  4  from  >  n^  -  2.0  to  n^  >  n^  -  4.5.  The  relative 
order  of  magnitude  and  the  general  sinusoidal  shape  is  similar 
to  the  drag  coefficient  for  the  SUBOFF  model.  Thus,  as  the 
body  becomes  shorter  and  rounder,  the  drag  force  increases. 
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The  Wave  Making  Lift  Coefficient  versus  Froude  Number  plot 
(Figure  23)  also  shows  an  oscillatory  character.  Two  local 
maxima  are  evident,  with  the  second  peak  being  of  greater 
magnitude  than  the  first.  As  the  exponents  are  raised,  the 
lift  coefficient  is  increased.  Also,  the  speed  at  the  peak  is 
shifted  slightly  to  the  right.  Of  special  note  is  the 
decrease  in  magnitude  above  F  -  0.43,  with  a  change  in  sign 
above  F  «  0.56.  As  seen  in  the  SUBOFF  model,  increasing  speed 
above  this  point  results  in  a  downward,  vice  an  upward,  sea- 
suction  force.  Overall,  the  shorter,  rounder  hulls  experience 
a  greater  lift  force  (almost  double)  than  do  the  longer, 
thinner  hulls. 

Likewise,  the  Wave  Making  Moment  Coefficient  versus  Froude 
Number  plot  (Figure  24)  is  oscillatory  in  nature.  With  the 
exception  of  a  few  low  speed  points,  all  values  are  negative, 
indicating  a  bow-down  or  stem-up  moment.  The  magnitude  of 
the  coefficient  increases  as  speed  is  raised  to  about  F  -  0.5 
indicating  that  a  higher  moment  is  produced  at  higher  speeds. 
But  e±»ove  F  -  0.5,  the  magnitude  slowly  drops.  The  shorter, 
rounder  hulls  produce  a  greater  bow-down  moment  than  do  the 
longer,  thinner  hulls.  Thus,  it  is  easier  to  keep  the  stem 
submerged  with  a  longer,  thinner  model. 

Figures  25  to  27  detail  the  incident  wave  mns.  The 
Incident  Wave  Drag  Coefficient  versus  Nondlmensional 
Wavelength  plot  (Figure  25)  shows  that  the  magnitude  decreases 
as  the  wavelength  increases.  Doubling  the  wavelength  drops 
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the  drag  coefficient  by  a  factor  of  two.  Thus,  as  the  waves 
become  longer,  the  drag  on  the  body  is  reduced.  For  the  two 
lowest  exponents,  the  drag  is  a  positive  value  (opposing 
force)  while  the  larger  esqponents  are  negative  (aiding  force) . 
The  more  streamlined  bodies  (lower  esqponents)  will  experience 
resistance  to  their  motion  while  the  shorter,  roiinder  bodies 
will  actually  be  pushed  in  the  longitudinal  direction. 

Figure  26  details  the  Incident  Wave  Lift  Coefficient 
versus  Nondimensional  Wavelength  plot.  The  lift  coefficient 
is  positive  throughout  (sea-suction  force) ,  peaks  at  a 
wavelength  of  ed)Out  one,  then  linearly  decreases.  As  the 
exponent  is  increased,  the  lift  coefficient  rises  through  all 
wavelengths.  However,  the  change  is  relatively  small  as  the 
lift  coefficient  varies  only  cd^out  25  percent  from  the  lowest 
to  highest  esqponent.  Thus,  the  shorter,  rounder  bodies 
e3q)erience  a  somewhat  greater  lift  force. 

Finally,  the  Incident  Wave  Moment  Coefficient  versus 
Nondimensional  Wavelength  plot  (Figure  27)  displays  a  similar 
shape  to  the  lift  coefficient  plot.  The  moment  coefficients 
are  positive,  indicating  a  bow-up  moment.  As  the  exponent  is 
raised,  the  moment  coefficient  is  increased.  About  a  20 
percent  increase  is  noted  from  the  lowest  to  the  highest 
exponent.  The  four  highest  coefficients  form  a  common  curve 
eibove  a  wavelength  of  1.5.  Thus,  the  shorter,  rounder  body 
will  feel  a  greater  bow-up  moment. 
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Varying  the  length  while  maintaining  the  displacement  and 
maximum  dicuneter  constant  result  in  changes  to  each  of  the 
three  coefficients.  Combining  the  effects  of  the  wave  making 
and  incident  wave  runs,  the  following  observations  are  made 
for  the  shorter,  roxinder  bodies  (higher  exponents) .  They 
produce  greater  drag,  which  is  due  mainly  to  the  wave  making 
effects.  The  Incident  waves  make  these  bodies  more 
susceptible  to  sea- suction  forces.  Also,  the  largest  bow-down 
moments  are  presented  to  these  bodies.  The  second  parametric 
study  examines  the  role  of  a  varying  dicuneter. 
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PARAiVlETP.IC:  DRAG  vs  FROUDE  NUMBER  (VARIABLE  LENGTH) 
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FRCUCE  NUMBER 


Figure  22:  Parametric  Study  of  Wave  Making;  Drag 
Coefficient  versus  Froude  Number  (Varicdjle 
Length) 


LIFT  COEFFICIENT 


^  q-3  PARAMETP.IC:  UFT  vs  rnCUDE  NUMBER  (VARIABLE  LENGTH) 


Figure  23:  Parametric  Study  of  Wave  Making:  Lift 

Coefficient  versus  Froude  Number  (Variable 
Length) 
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X  iq-^PARAMETPIC:  MOMENT  vs  FP.CUOE  NUMBER  C^ARIABLH  LENGTH) 


Figure  24:  Parametric  Study  of  Wave  Maiking:  Moment 
Coefficient  versus  Froude  Number  (Variable 
Length) 
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PAP^METP.IC  INC  WAVE  UFT  vs  LAMBDA  (VARIABLE  LENGTH) 


Figure  26:  Parametric  Study  of  Incident  Waves;  Lift 
Coefficient  versus  Nondimens ional  Wavelength 
(Variable  Length) 
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Figure  27:  Parametric  Study  of  Incident  Waves:  Moment 
Coefficient  versus  Nondimens ional  Wavelength 
(Variable  Length) 
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C.  EFFECTS  OF  DIAMETER 


In  the  second  parametric  study,  the  displacement  again 
remains  constant  at  24.692  ft^,  emd  the  overall  length  and 
depth  are  fixed  at  14.3  ft  and  1.43  ft,  respectively.  In  each 
run,  the  diameter  is  the  variable  and  the  effects  on  the  drag, 
lift  and  bow-up  moments  are  examined.  The  diameter  for  each 
is  obtained  by  finding  the  roots  of  (90)  for  D.  Again  the 
values  of  es^onents  determine  the  prismatic  coefficients  (see 
Tcible  II).  The  root  closest  to  1.667  (the  original  diameter) 
is  chosen  for  the  body  diameter.  Table  III  presents  the 
exponents  and  their  corresponding  diameters. 


TABLE  III.  EXPONENTS  AND  CORRESPONDING  DIAMETERS 


Ha  -  n. 

DIAMETER  (ft)  I 

2.0 

1.78  1 

2.5 

H 

• 

H 

3.0 

1.66 

3.5 

1.63 

4.0 

— nr— 1 

L_ll_ 

...  J 
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As  the  Table  III  clearly  shows,  the  larger  exponents  result  In 
smaller  diameters  to  maintain  the  same  displacement.  As  the 
ends  become  more  rounded,  the  mcocimum  diameter  must  be 
reduced.  The  results  generated  by  the  con^uter  program  are 
shown  in  Figures  28  to  30  for  the  wave  making  runs  and  in 
Figures  31  to  33  for  the  incident  wave  mns. 

An  oscillatory  set  of  curves  is  obtained  for  the  Wave 
Making  Drag  Coefficient  versus  Froude  Number  plot  (Figure  28)  . 
Two  peaks  are  observed,  with  the  peaks  shifting  to  higher 
speeds  as  the  exponents  increase.  As  the  exponents  are 
raised,  the  magnitude  of  the  drag  coefficient  increases,  and 
is  more  than  doubled  from  the  lowest  to  the  highest  exponent. 
As  the  body  becomes  rounder  on  the  ends  and  thinner  in  the 
middle,  the  drag  force  produced  on  the  body  increases. 

Figure  29  details  the  Wave  Making  Lift  Coefficient  versus 
Froude  Number  for  each  exponent.  The  oscillatory  nature  is 
again  apparent,  but  the  shapes  of  each  curve  shows  larger 
variations  at  lower  speeds.  For  increasing  exponents  the 
first  peak  shifts  slightly  to  the  right  and  the  magnitude 
increases  while  the  second  peak  shifts  to  higher  speeds  and 
the  magnitude  decreases.  As  in  the  varying  length  study, 
cibove  F  0.55,  all  lift  coefficients  drop  below  zero.  At 
lower  speeds,  the  body  produces  an  upward  force  and  at  higher 
speeds  a  downward  force  is  produced.  The  rounder,  thinner 
(greater  exponents)  body  produces  the  greatest  lift  only  in  a 
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small  region  about  F  «  0.3.  In  the  other  lower  speed  regions 
(less  than  F  -  0.55),  the  less  rounded,  fatter  (lower 
esqponents)  body  produces  the  higher  sea- suction  force.  Thus, 
to  reduce  the  effects  of  sea- suction  forces,  a  rounder, 
thinner  body  would  be  prefered^le.  Once  again,  higher  speeds 
mitigate  the  sea- suction  forces. 

The  last  plot  (Figure  30)  shows  the  Wave  Making  Moment 
Coefficient  versus  Froude  Nximber  for  various  e3q>onents.  Once 
again,  the  oscillatory  nature  and  the  negative  (bow-down) 
moment  coefficients  are  present.  A  tighter  set  of  curves  is 
obtained  than  for  the  varied)le  length  case.  At  lower  speeds 
(less  than  F  >  0.35),  the  higher  exponents  yield  the  greatest 
coefficient  magnitudes,  while  the  higher  speeds  (F  »  0.35  to 
F  «  0.6)  result  in  lower  coefficient  magnitudes.  The  less 
rounded,  fatter  bodies  result  in  almost  no  moment  at  low 
speeds,  but  produce  a  larger  moment  as  speed  was  increased. 
The  rounded,  thinner  bodies  produce  a  more  uniform  moment 
throughout  the  range  of  speeds  encountered.  Overall,  every 
body  form  experiences  a  bow- down  moment. 

Figures  31  to  33  detail  the  Incident  wave  runs  as  the 
diameter  is  varied.  The  Incident  Wave  Drag  Coefficient  versus 
Nondimens ional  Wavelength  plot  (Figure  31)  shows  a  decreasing 
magnitude  coefficient  throughout  all  wavelengths.  Similar  to 
the  varying  length  case,  the  the  lowest  exponents  yield 
positive  coefficients,  while  the  highest  exponents  produce 
negative  drag  coefficients.  As  the  body  becomes  rounded  and 


91 


narrow  (higher  exponent) ,  the  incident  waves  actually  aid  the 
body's  motion  instead  of  opposing  it.  Still,  the  effects  are 
very  small  when  compared  to  the  wave  making  case. 

Figure  32  details  the  Incident  Wave  Lift  Coefficient 
versus  Nondimens ional  Wavelength  plot.  A  peak  near  a 
wavelength  of  one  is  seen  with  a  steady  decline  thereafter. 
As  the  exponents  increase,  theposition  lift  decreases  slightly 
(odDout  five  percent) .  This  shows  that  a  thinner,  rounder  body 
will  experience  a  smaller  sea- suction  force  than  a  wider, 
streamlined  body. 

The  last  plot.  Figure  33,  displays  the  Incident  Wave 
Moment  Coefficient  versus  Nondimens ional  Wavelength.  A 
similar  shape  to  the  varying  length  plot  is  observed.  As  the 
exponents  increase  from  lowest  to  highest,  the  positive  (bow- 
up)  moment  decreases  by  about  40  percent.  As  with  the  lift 
coefficient,  the  thinner,  rounder  body  will  feel  a  smaller 
bow- up  moment  than  a  wider,  streamlined  body. 

In  summary,  changing  the  maximum  dicuneter  and  the  body 
shape  for  a  given  length  and  displacement  modify  each  of  the 
three  coefficients.  The  less  wider,  streamlined  (lower 
exponent)  body  produces  the  lowest  drag,  the  greatest  sea- 
suction  force,  and  the  most  variadjle  bow- down  moment. 
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PARAMETRIC:  DRAG  vs  FROUCE  NUMBER.  (VARIABLE  DIAMETER) 
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Figure  29:  Parametric  Study  of  Wave  Making:  Lift 

Coefficient  versus  Froude  Number  (Variable 
Diameter) 
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Figure  30:  Parametric  Study  of  Wave  Making: 

Coefficient  versus  Froude  Number 
Diameter) 
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Figure  31:  Parametric  Study  of  Incident  Waves:  Drag 

Coefficient  versus  Nondimens ional  Wavelength 
(Variable  Diameter) 
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Figure  32:  Parametric  Study  of  Incident  Waves;  Lift 
Coefficient  versus  Nondimens ional  Wavelength 
(Variable  Diameter) 
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Figure  33:  Parametric  Study  of  Incident  Waves:  Moment 
Coefficient  versus  Nondimens ional  Wavelength 
(Variable  Diameter) 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  work  examined  two  significant  free  surface  effects  in 
a  submerged  body.  The  waves  produced  by  the  moving  body 
interact  with  the  free  surface,  euid  in  turn,  produce  forces 
euid  a  moment  on  the  body.  Secondly,  the  incident  waves  on  the 
otherwise  calm,  free  surface  also  produce  its  own  forces  and 
moment  on  the  body.  The  combination  of  these  two  effects 
produce  a  total  response  on  the  body.  The  wave  making  effects 
were  the  major  contributor  to  the  body  drag,  causing  a  force 
which  resisted  the  body's  longitudinal  motion.  The  incident 
waves  produced  the  largest  upwards  lift  force  on  the  body, 
pulling  the  body  to  the  free  surface.  A  combination  of  the 
wave  making  and  incident  wave  effects  were  responsible  for  the 
total  bow- down  moment  on  the  body.  Understanding  these 
effects  will  aid  the  designer  in  optimizing  future  submarine 
designs  and  will  assist  submarine  crews  in  effectively 
maintaining  depth  control  at  periscope  depth.  On  the  design 
side,  the  data  from  this  work  could  be  used  directly  in  the 
simulation  design  based  process  or  in  design  of  an  automatic 
depth  control  system.  On  the  operational  side,  the  data  could 
be  input  into  training  simulators  to  better  model  the 
submarine  operating  at  periscope  depth.  Additionally,  the 
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results  could  be  incorporated  into  the  existing  Ship's  Fleet 
Mission  Program  Library  (SFMPL)  to  aid  the  crew  in  trimming 
the  submarine  prior  to  proceeding  to  periscope  depth. 

To  further  inprove  the  knowledge  of  free  surface  effects 
on  submerged  bodies,  the  following  recommendations  are 
proposed.  First,  the  viscous  effects  in  the  body's  boundary 
layer  can  be  modelled  to  yield  additional,  significant  forces 
etnd  moments.  These  effects,  especially  flow  separation  from 
the  body  and  its  appendages  would  greatly  add  to  the  potential 
flow  results  discussed  in  this  work.  Secondly,  the  effects 
from  the  bottom  or  from  other  nearby  sources  (proximity 
effects)  would  likewise  increase  the  accuracy  of  the  resulting 
forces  and  moments.  Finally,  the  effects  of  incident  wave 
diffraction  (which  was  neglected  in  this  work)  would  produce 
additional  forces  and  moments.  Also,  modelling  the  seaway  as 
a  random  wave  vice  a  sinusoidal  wave  would  produce  a  more 
realistic  seawave.  These  recommendations  would  serve  to 
enhance  the  results  presented  in  this  work. 
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